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1. Introduction 

A substantial increase in the ratio of strange over non-strange particles between 
proton-proton and central nucleus-nucleus collisions has been observed experimen¬ 
tally both at CERN and AGS [1]. Such an increase was expected as a consequence 
of quark gluon plasma (QGP) formation due to an increase of the strange quark 
production rate in the deconfined phase [2]. However, as in the case of other possible 
signals of QGP, one has to check whether such a phenomenon is typical of a central 
nucleus-nucleus collision or, on the contrary, it is already present in proton-nucleus 
interactions *. Unfortunately, the experimental situation is not conclusive. While 
in the data from the NA35 collaboration strangeness is enhanced between pp and 
pS, the effect is not seen when taking an average over three different nuclear targets 
[4]. However, two other experimental facts indicate that strangeness enhancement 
is not only present in central nucleus-nucleus collisions. Indeed in pp collisions, the 
ratio K/tt increases both with increasing energy and with increasing multiplicity [5]. 
Also the ratio A/p in pp collisions increases from 0.26 at ^/s = 20 GeV [6] to 0.39 
at ^/s = 1 800 GeV [7]. Furthermore, results from the NA36 collaboration [8] show 
a large strangeness enhancement between peripheral and average nucleus-nucleus 
collisions. 

Actually, it has been emphasized [9-11] that strangeness enhancement is ex¬ 
pected to occur in independent string models such as the Dual Parton Model (DPM) 
[9] and the Quark Gluon String Model (QGSM) [12]. Indeed, the average number 

* For instance, in the case of J/i/ suppression, R. Salmeron, followed by many 
other authors, has shown that the effect could be explained by extrapolation of the 
corresponding hadron-nucleus data [3]. 
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of collisions per participant nucleon increases between pp and pA as well as be¬ 
tween pA and central AB collisions. The resulting extra particles are produced in 
the fragmentation of strings involving sea constituents at their ends. Since strange 
sea quarks are present in the nucleon sea, the ratio of strange over non-strange 
secondaries will also increase. The value of this increase depends on the ratio 

^ ^ 2(g + g) 

u + d + u + d 

in the nucleon sea. We shall use a value of S in the range 0.4 to 0.5, consistent 
both with strangeness suppression in string fragmentation at the highest available 
energies [7] and with conventional parton distributions for DIS [13]. In order to 
explain the observed enhancement of A and A, two approaches have been proposed. 
In DPM, diquark-antidiquark pairs have been introduced in the nucleon sea - with 
the same ratio, relative to quark-antiquark pairs, as in the string breaking process 
[10, 11]. The second approach is string fusion [14, 15]. However, both mechanisms 
have a common caveat, namely, they lead to the same increase in absolute value of 
strange baryons and antibaryons - while experimentally the excess of A’s is much 
larger than that of A’s. 

In the present work we solve this problem by introducing final state interactions 
of co-moving pious and nucleons which increase the number of produced A’s via 
the binary interactions tt + N K + A {K* -|- A, ...). This is a fundamental 
mechanism in hadron gas models [2] and has been incorporated into some Monte 
Carlo codes [12, 15, 16]. Clearly, such an interaction is very efficient in producing 
A’s - since the density of both pions and nucleons is high as compared to that of 
A’s. Actually, it was shown in Ref. [11], that using the experimentally measured 
pion and nucleon densities in central SS collisions, together with the cross-section 
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of the above reaction (which is about 1.5 mb at the maximum near threshold), 
one precisely obtains the missing number of A’s. Obviously, the conjugate reaction 
TT + N ^ K + A {K* + A) will also produce some extra A but the absolute yield is 
much smaller than the one of A’s due to the smallness of the N density as compared 
to that of A^’s. 

In this paper we compute, in the framework described above, the rapidity 
distributions of A and A in SS and Pb-Pb collisions. A substantial increase of the 
ratios A/'k~ and A/tt~ (70 % and 50 %, respectively) is predicted between SS and 
Pb-Pb. 

Extra evidence for the necessity of a final state interaction is obtained from 
baryon stopping. The rapidity distribution of proton minus antiproton in central 
SS collisions, computed in independent string models, has a dip at y* ~ 0 which 
is much more pronounced that the one observed experimentally [17] - while the 
corresponding distributions both in pp and in peripheral SS collisions are well re¬ 
produced. In other words, the stopping power in the model is too small as compared 
to experiment and rescattering of the produced baryon is required in order to in¬ 
crease it [12] [15]. 

The plan of the paper is as follows : in Section 2 we give the general formulae 
of DPM for nucleus-nucleus collisions in the presence of diquark-antidiquark pairs 
in the nucleon sea and compute the rapidity distributions of A and A in pp and SS 
collisions without final state interactions. In Section 3, we introduce the final state 
interactions and give the modified A and A rapidity distributions in S'S' collisions 
as well as the model predictions for Pb-Pb collisions. The rapidity distribution 
of pious and protons in pp, SS and Pb-Pb collisions are also given. A list of 
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average multiplicities for pious, kaous, baryous aud autibaryous without aud with 
fiual state iuteractious is giveu iu Table 1. Sectiou 4 coutaius some fiual remarks 
aud couclusious. Iu Appeudix 1 we give the momeutum distributiou fuuctious aud 
fragmeutatiou fuuctious used iu this work, as well as the formulae for the hadrouic 
spectra of the iudividual striugs. 

2. The Model 

Iu DPM aud QGSM, the domiuaut coutributiou to multiparticle productiou iu 
pp collisious is a two striug productiou - the two striugs beiug stretched betweeu 
the valeuce diquark of a protou aud a valeuce quark of the other protou, aud vice- 
versa. Iu fact, iu this paper we shall restrict ourselves to this compoueut (see 
Sectiou 4 for a discussiou ou this approximatiou). Therefore we have exactly two 
striugs per iuelastic collisiou. Wheu the uumber of iuelastic collisious experieuced 
by a uucleou iu a uuclear collisiou is larger thau oue, the extra striugs (two per 
each extra collisiou) have sea quark coustitueuts at the eud correspoudiug to the 
multiply wouuded uucleou. As explaiued iu the lutroductiou, siuce strauge sea 
quarks are preseut iu the uucleou sea, the ratio of strauge over uou-strauge particles 
will iucrease. Oue cau object, however, that if the fractiou S of strauge quarks 
(defiued iu the lutroductiou) is the same as iu the striug breakiug process, theu the 
striugs with sea quarks at their euds will produce the same strauge over uou-strauge 
ratio as iu the striug breakiug process - with uo uet straugeuess euhaucemeut. While 
true at asymptotic euergies, this is uot the case at preseut oues. For iustauce, let us 
cousider the K~/7t~ ratio. Iu NN collisious at 200 GeV/c, the experimeutal value 
of this ratio is about 0.05. Siuce the values of S we cousider are much larger, it is 
clear that the ratio K~Itx~ will iucrease. Auother example is the ratio A/p which. 
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as discussed in the Introduction, reaches a value close to 0.4 at Fermilab [7]. In the 
following we shall use the value S = 0.5. However, taking S = 0.4 changes very 
little our results [11]. 


The general formulae for the single particle inclusive rapidity distribution in 
nucleus-nucleus collisions is the following [9] : 

1 


dy 


O'AB 


E 


cr, 


nA,nB,n 


ub - ua) riA ( 


nA[NZ^-f{y) + Nf-^'^" 


IV 'rt'rlv 


(»)) + 


nA,nB,n 


(ns - ha) (nZJb (y) + nZJZ ) + (^ - nB) (n^Jb iv) + ^ZJb (y) 


+sym {ha ^ ^b)} 


( 2 . 1 ) 


Here ^ is the cross-section for n inelastic nucleon-nucleon collisions involving 

Ha and ns participating nucleons, gab = On,nA,nB inelastic nucleus- 

nA,nB,n 

nucleus cross-section, and fiA = n/uA, fJ^B = n/nB- N{y) are the inclusive spectra 
of each individual string. Note that the total number of strings in Eq. (2.1) is 2n 
(two for each nucleon-nucleon collision). 


In the following we restrict the discussion to the case A = B. Eq. (2.1) can be 
written in a simple, albeit approximate, form 

dNApAT 


dy 


^nANrAlAy) + N£7f(y) 


+(n - nA ( n£;A (y) + 


( 2 . 2 ) 


Eq. (2.2) is obtained by neglecting the dependence of the string densities N 
on the /u’s (which are replaced by their average values /f). The summations in (2.2) 
are then trivial - the numbers ua and n being replaced by their average values 
fiA and n. The approximate Eq. (2.2) is quite good at mid-rapidities. At large y 
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it becomes worse. However, in the rapidity region where the baryon densities are 
sizeable (y* < 2.5), the errors involved are 5 to 10 %. 

In Eqs. (2.1) and (2.2) only Qs-Qs pairs in the nucleon sea are considered. As 
discussed in the Introduction, we now assume that qqs-QQs pairs are also present with 
relative fraction a. The modification in Eqs. (2.1) and (2.2) are straightforward. 
Eq. (2.2) is changed into : 


dNApAT 

dy 




Ma ,/J-A 


(y)) + (n - ua) I (1 - 2a) ((y) 


+N£/:(y)) + a (N<i,f(y) + (v) 


(2.3) 


In the following, we compute the rapidity distributions of A and A in nucleus- 
nucleus collisions using Eq. (2.3). The rapidity distributions of the individual 
strings are obtained by convoluting momentum distribution functions and fragmen¬ 
tation functions. All relevant formulae (including absolute normalization) are given 
in Appendix 1. 


In numerical calculations we use the value a = 0.1 taken from JETSET [18]. (A 
somewhat smaller value, a = 0.07, is obtained in ERITIOE [16] and in BAMJET 
[19]). The values of ua and n are obtained in the optical approximation of the 
Glauber model, with Saxon-Woods profiles, for impact parameter \b\ < 1 fm. We 
get in SS collisions : ua = 31, h = 55 and in Pb-Pb collisions ua = 192, n = 676. 


Before giving the numerical results, let us make some remarks concerning 
Eq. (2.3). Its first two terms are the same as in NN collisions - except for the 
value of flA - which is 1 for NN, 2 for central SS and 4.5 for central Pb-Pb. The 
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third and fourth terms give negligeably small contributions to baryon-antibaryon 
production due to the large, {tua + my)^, thresholds in the short Qs-Qs strings. The 
remaining terms contain the contribution of sea diquarks and produce A and A by 
pairs. Detailed formulae for the rapidity distributions of baryons and antibaryons 
resulting from the various terms of Eq. (2.3) are given in Appendix 1. 

In Fig. 1 we give the Feynman x distributions of A and A in pp collisions. 
They are obtained from the first two terms of Eq. (2.3) (i.e. from the sum of 
the contributions of the two qq-q valence-valence strings) with = 1. As shown 
previously in Ref. [20], the agreement with available data is quite good. Note that 
the shape of the fragmentation functions is determined from Regge intercepts (see 
Ref. [20] and Appendix 1). 

The corresponding rapidity distributions in central SS collisions are given in 
Figs. 2 and 3. We see that the yields of both A and A are smaller than the data. The 
A yield, in particular, is considerably smaller than the experimental one. There¬ 
fore, we need a mechanism to produce extra A’s. As already mentioned in the 
Introduction this can be achieved by a final state interaction of co-moving secon¬ 
daries 7T + N ^ K + A {K* -|- A, ...). A smaller number of A will also be produced 
via TT + N ^ K + A {K* + A, ...). In the next section we discuss in detail these 
final state interactions. 

3. The Final State Interaction of Secondaries 

In independent string models, it is assumed that secondaries produced in dif¬ 
ferent strings are independent from each other. Introducing a final state interaction 
of secondaries represents a departure from the independent string picture. In fact. 
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it corresponds to the simplest form of string interaction. It should be stressed that 
it is not known how to treat quantitatively such an interaction. Here we adopt 
the formalism introduced in Ref. [21]. Note that, at moderate energies, such a 
treatment of string interaction can be justified in the framework of Reggeon Field 
Theory [22]. 

Let us consider the number of A’s produced by unit of space time volume d^x 
resulting from the interaction tt + > RT + A {K* + A, ...). One has [21] 


dN^ 


= < (7 > pT,{x) pp{x) 


(3.1) 


where p{x) are space-time dependent particle densities and < cr > is the interaction 
cross-section - properly averaged over the momentum distributions of the colliding 
particles. We take < a >= 1.5 mb which is the value oi tt~ + p ^ + A/S° at 

its maximum near threshold. Beyond the threshold these two-body cross-sections 
decrease rapidly but quasi two-body interactions convert this sharp decrease into a 
mild increase. 


We use cylindrical space-time variables : the longitudinal proper time is r = 


t-\-z 

t — Z 


. In these variables d^x = 


\/t‘^ — z‘^ and the space-time rapidity y = 

T dr dy d^s, where s is the transverse coordinate and d‘^s an element of transverse 
area. As customary, we assume boost invariance - i.e. that the space-densities Pi{x) 
are independent of y - and assume, furthermore, that the dilution in time of the 
density is mainly due to longitudinal motion. With this two assumptions we have 
[ 21 ] : 


Pi{r,y,s) = piir.s) — 

T 


(3.2) 
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Using (3.2) it is trivial to integrate Eq. (3.1) on the variable r. One gets 


dNA 

dy 


= / d'^s 


,2 dN-u: dN. 


Ujiy qr UjI y n ^ r / \ / T 

, ,2 ^ ^2 3 < (7 > £n (r + ro)/ro 

dy d^s dy d^s 


(3.3) 


where tq is the formation proper time and r the time during which the final in¬ 
teraction takes place. We take tq = 1 fm and r = 3 fm - based on interferometry 
measurements which indicate a very short time of particle emission and a freeze out 
time To + r ~ 4 fm [23]. The factor 3 results from the product of three pion times 
two nucleon species, divided by a factor 2 - which is due to the fact that one out of 
two 7tN combinations have a negligeably small cross-section. For instance among 
the TT^p combinations only n~p has to be considered. 

Let us now explain how to compute dNi/dy d?s in Eq. (3.3). The formulae 

given in Section 2 (see for instance Eq. (2.3)) give the expression of dN/dy - which is 

the result of integrating dN/dy d?s over d?s. Actually, the only dependence on s is 

in the average values ha, ns and n. In the Glauber model, these average quantities 

—^ 

are known as a function of b and s. In the optical approximation one has 

n(b, s) = ^TA(b-g) Tb{s) , (3.4) 

O'AB 


and 



A 

O'AB 


Ta(6 - ^ o-jvAis) 


where 

0‘N/l(s) = 1 - (1 - (T 


(3.5) 


Here a is the nucleon-nucleon inelastic cross-section and TA{b) the nuclear profile - 
for which we use the standard Saxon-Woods expression. Clearly, when computing 
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the r.h.s. of Eq. (3.3) one has to perform the products n\(b,^, nA{b,^ 

n{b, ^ and integrate them over both (Ps and d?b. It is easy to see that such integrals 
of products are smaller than the product of integrals dNT^/dy times dNp/dy and 
that the ratio of the former over the latter is roughly proportional to ~ . 

Therefore, the number of A’s resulting from the final state interaction, will increase 
quite fast with increasing A. Indeed, at mid-rapidities, the product dNj^/dy times 
dNp/dy increases with A approximately as A^^^ and thus the number of A’s increases 
approximately as A^. 

In order to compute the rapidity distributions of A resulting from the final 
state interaction (3.1) one has to compute first the densities of pions and protons. 
This has been done using the formulae in Appendix 1. 

The rapidity distributions of h~ in NN, SS and Pb-Pb collisions are given 
in Figs. 4 and 5. The agreement with experiment is reasonable. Note that, due 
to the two string approximation in NN, the rapidity distribution tends to be too 
broad. Indeed, the contribution of the Qs-Qs strings, although small at 200 GeV/c, 
is concentrated at mid-rapidities. This is partly compensated by choosing steeper 
diquark fragmentation functions (see Appendix 1). In any case, our aim here is 
not to get the best possible description of h~ distributions, but only to have a 
reasonable form of dN^-/dy in order to use it in Eq. (3.1). Following Ref. 4, we 
have taken dN^-/dy = 0.93 dN^- /dy. 

The Feynman x distributions of proton and antiprotons in pp collisions are given 
in Fig. 6. The agreement with data is quite good. Note that diffractive dissociation, 
which produces a peak near x = 1 in the p distribution has not been included in 
the model. The corresponding results for peripheral and central SS collisions are 
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given in Fig. 7. Actually, both the data and the theoretical curves correspond to 
the difference p — p. The shape of the rapidity distribution is well reproduced in the 
model for peripheral collisions but not for central ones. In the latter case, the shape 
of the experimental distribution is much flatter than the theoretical one. Note, 
however, that the average multiplicity is the same for the two distributions. This 
disagreement between an independent string model prediction and experiment is 
very significant and we are going to discuss it in some detail. First of all, in DPM 
the stopping power is entirely controlled by the momentum distribution functions. 
More precisely, when the average number of collisions per nucleon increases (i.e. 
when pA in Eq. (2.3) increases), the maximum of the proton spectrum is shifted to 
smaller rapidities (see Fig. 2 of Ref. [24]). This effect is seen in Fig. 7. However, 
in central SS collisions pA and, therefore, the effect is small. Moreover, the 
decrease of the proton spectrum between its maximum and y* ~ 0 is very steep even 
for a large number of inelastic collisions. Thus, the dramatic flattening of the proton 
rapidity distribution between peripheral and central SS collisions in Fig. 7 cannot 
be reproduced in the model, and, as discussed in the Introduction, rescattering of 
the produced nucleons is needed. 

We can now compute the rapidity distribution of A resulting from the final 
state interaction tt + N ^ K + A {K* + A, ...). The results are shown in Fig. 2. 
The agreement with experiment is quite good. 

The corresponding results for A production, including the effect of the final 
state interaction tt + N ^ K + A {K* + A, ...), are given in Fig. 3 (dashed line). 
We see that the final state interaction has produced a substantial increase of the 
A yield, which is now in qualitative agreement with experiment. We conclude that 
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the main features of A and A production in central SS collisions can be understood 
in the scenario presented in this paper. 

We turn next to central Pb-Pb collisions at 160 GeV/c per nucleon. In this 
case the rate of strange baryon production from the final state interactions is 
very high at mid-rapidities due to the large value of the pion density and one 
can not ignore any longer the crossed reactions tt -|- A —> K + N and tt -|- A — 
K + N. Indeed, it turns out that, at mid-rapidities, equilibrium between strange 
and non strange baryons is reached before the freeze out time tq + r (assumed to 
be the same as in central SS collisions). When this equilibrium is reached one has 
< n >p + < n > n = 1-6 < n >a (where 0.6 < n >a takes into account the 
baryons) and the strange baryon multiplicity does not increase any longer. The 
predictions for the rapidity distribution of A and A in central Pb-Pb are shown in 
Figs. 2 and 3 (upper lines scaled down by a factor 5). 

A close look at Fig. 3 shows that the computed A distribution in SS collisions 
has a shallow dip at y* = 0 which is not seen in the data. Probably this is related to 
the lack of baryon stopping discussed above and requires further rescattering of the 
produced A’s. In order to have a rough estimate of the effect of the baryon rescat¬ 
tering, we have introduced some extra stopping in the proton spectrum (without 
changing the proton average multiplicity) by shifting it by Ay = 0.5. The modified 
p — p rapidity distribution is shown in Fig. 7 (full upper line). Using such shifted 
proton distributions in Eq. (3.3), we obtain a flatter A distribution both in SS and 
Pb-Pb collisions (full lines in Fig. 3), which for SS collisions is in good agreement 
with experiment. 


Finally, the average multiplicities of pions, kaons, protons lambdas and anti 
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lambdas in SS and Pb-Pb collisions are given in Table 1. (A discussion of cascades 
and omega production in the present framework can be found in Refs. 10, 11). We 
observe that the ratio K/h~ {A/h~) increase by about 50 % (70 %) between SS 
and Pb-Pb collisions. 

4. Conclusions 

We have studied meson and baryon strangeness enhancement in an independent 
string model with quark-antiquark and diquark-antidiquark pairs in the nucleon sea. 
The expe rimental data have forced us to depart from string independence and to 
introduce a final state interaction of co-moving secondaries. The necessity of final 
state rescattering also follows from the study of the shape of the proton rapidity 
distribution in central SS collisions. Final state interactions are the simplest forms 
of string interaction - which is a first step in order to drive the system toward 
equilibrium. Using all these ingredients, we have obtained a reasonable agreement 
with experiment. We feel that a more accurate description is not possible in the 
present state of the art. 

The present model has many uncertainties that we want to discuss briefly. 
Apart from the poor knowledge of the parameter S and a, there is an uncertainty 
resulting from the hadronization of short strings and, in particular, of strings in¬ 
volving sea diquarks (see Appendix 1). Furthermore, the introduction of final state 
interactions rises a lot of questions. Not only a rigourous treatment of the final 
state interaction does not exist, but, moreover, we can not be sure that other re¬ 
actions (for most of which the cross-sections are not even known) do not affect 
our results. (In particular we have neglected annihilation reactions of nucleons and 
antilambdas or antinucleons and lambdas which have a very large cross-section at 
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threshold. However, this cross-section drops very fast after the threshold and the 
corresponding value of < cr > in Eqs. (3.1), (3.3) is not known). Despite all these 
drawbacks, the effects discussed in the present work are important and cannot be 
ignored in strangeness enhancement studies. 

A characteristic feature of the model presented here is that it predicts a contin¬ 
uous strangeness enhancement, not only from proton-proton to proton-nucleus and 
to average (i.e. minimum-bias) nucleus-nucleus collisions, but also in proton-proton 
when energy or event multiplicity is increased. Of course to get this latter effect it is 
necessary to give up the two string approximation in pp collisions. This assumption 
is of no consequence for the results presented here, provided a good description of 
pp data is achieved [25]. Indeed the results for nucleus-nucleus collisions are deter¬ 
mined by the input in pp (or rather NN) collisions. However, when multistrings 
are introduced, one obtains [26] an increase of the ratio K/tt in pp collision when 
increasing either the energy or the event multiplicity, as observed experimentally. 
Moreover, the increase with energy [7] of the ratio A/p, can also be understood in 
this way. Note, however, that the contribution to strangeness enhancement result¬ 
ing from the final state interaction (which is very important for A production) is 
most effective in central nucleus-nucleus collisions where the densities of secondaries 
are very large. 

In conclusion, the study of strangeness enhancement is a very active field. This 
is clear from the large number of recent work on the subject [1]. Actually, during 
the latest stages of the present work [27] many interesting papers on this subject 
have appeared in the literature [28-31]. Some of them give also predictions for Pb- 
Pb collisions. The theoretical ideas in these papers and their predictions are, in 
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general, quite different from ours. Forthcoming results with a lead beam at CERN 
will certainly trigger more activity and give new insight into this very interesting 
field. Predictions for RHIC and LHC energies are also needed. 
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Appendix 1 


Hadron Spectra of Individual String 

Both in DPM and QGSM the hadron spectra of the individual strings are given 
by a convolution of structure functions and fragmentation functions. In DPM, one 
has for the rapidity density s) of hadron h in, say, a diquark-quark string : 

f dxi [ dX2 pIUxi) pI(x2) 


M1IM2 

dy 


'0 


'0 


(y 'Sstr.) di^Sgtr. ■Sfft.r.) 


(Al) 


Here and momentum distribution functions of diquark and quark 

which depend on the light cone momentum fractions Xi and X 2 - The index pi 
(P 2 ) denotes the number of inelastic collisions suffered by the nucleon to which the 
diquark (quark) belongs. The squared invariant mass of the string is Sgtr = X 1 X 2 S 
and A = 1/2 log[xi/x 2 ] is the rapidity distance between the C. of M. of the string 
and the overall C. of M. At fixed values of xi and X 2 the string hadronization is 
just given in terms of qq and q fragmentation functions : 


I for y > A 


dy 


■ (y ^5 ■Sstr. ) 


G'J(y) for y < A 


(A.2) 


where G{y) are fragmentation functions (see below). Of course /dy = 0 

for y outside the rapidity interval determined by the string endpoints. Note the 
presence of a threshold d-function in Eq. (A.2). The value of Sthr corresponds to 
the squared mass of the lightest system that can be produced in association with 
the triggered hadron h*. 


For instance for a qqv-qv string fragmenting into a A we have Sthr = (^t-a + 
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An important feature of DPM and QGSM is the fact that the structure function 
is determined in terms of Regge intercepts. In DPM, the structure function in the 
presence of 2/c such constituents is given by [9] : 

pk{xi ■ ■ ■ X 2 k) = C2k ^ ■ ■ ■ x~l_^xy ^5 ^ ■ (^-3) 

The index 1 refers to the valence quark and the index 2k to the diquark. All other 
indices refer to sea quarks and antiquarks. The constant C2k is obtained from the 
normalization of p2k to unity. The structure function of each individual constituent, 
which appears in Eq. (A.l), is obtained by integrating (A.3) over the variables of 
all other constituents. (In actual calculations the singular factor x~^ is replaced by 
x~^ with X = and p = 0.1 GeV ; varying p between 0.05 and 0.3 

produces only minor changes in our results). 

In QGSM the formulae are somewhat different [7]. The diquark and quark 
fragment independently from each other and the threshold resulting from the 0- 
function in (A.l) is not present. Moreover in Eq. (A.3) sea quarks are assumed 
to have, at present energies, a x~^l‘^ behaviour at x ~ 0 - the same as for valence 
quarks. In this case the integration of p in (A.3) over all variables except one can 
be performed analytically. One gets 

p‘g(x) = Cll - xp-^'^ . (A.i') 

ntiKY■ Eor that of a us-d string we have Sthr = Likewise for pion (kaon) 

production in a ud {sd) string we have Sthr = - where Pt^ is the pion 


transverse mass. 



19 


The coefficients C are determined (analytically) from the normalization to unity of 
p{x). Note that the x~^ behaviour in DPM strongly reduces the probability of a 
sea quark to cross to the opposite hemisphere (as compared to an x~^ behaviour, 
which corresponds to a fiat distribution in rapidity). As a consequence, the results 
obtained in DPM and QGSM are quite similar - at least at the energies considered 
here. 

In the first DPM papers, fragmentation functions were just derived from a fit to 
pp (or e'''e“) data. An important break through was achieved when it was realized 
that, not only momentum distribution functions, but also fragmentation functions 
could be derived in terms of Regge intercepts [7]. In this way the predictivity of 
the model was greatly increased. Indeed, it has been possible to predict the shape 
of the rapidity distributions of tt, AT, p, A, etc in pp collisions - with absolute rates 
as the only free parameters. 

In the following, we use the fragmentation functions of baryons and antibaryons 
obtained in this way. More precisely, the form of the fragmentation functions for A 
and A are given by [20] 

A = A = 1(1 - zf'’ 

Dt = A = 1(1 - zf’- 

A„,2 = Ai2 = Ad = 1(1 - z)"" 

A„,2 = Ad,2 = = (1 - . (.4.4) 

Here Dqq ^2 denotes the fragmentation of a diquark into a non-leading baryon. We 
have to add to it a (dominant) term, Dqq^i, corresponding to leading baryon pro- 
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duction, given by 


nA _ tjA 

^uuA — ^ 


ddA 


- Z)»-5 


(^.5) 


Turning to p and p production, we have [20] 


Dl = Dl= i(l - Zf 

DS = |(1-Z)" ; DS = A(1-Z)2 

0 L ,2 = DL .2 = - 0^,2 = ^(1 - Zf 

nP — nP — nP — 

^dd,2 ~ ^uu,2 — ^ud,2 ' 




(A6) 


and, for the leading proton. 


i'L.i = - Z)“ ; i?L,i = 4(l-^) 




(A.T) 


Finally for neutron production we have 




D'i = OS 


T^n _ T^p 
-‘-^uu,2 ^uu,2 ’ 


T)n _ 7-)P 

^udp{2) ^udp{2) 




(A.S) 


The values of all powers of Z and 1 — Z in Eqs. (A.4)-(A.8) are determined in 
terms of Regge intercepts*. For example, the power 0 of 1 — .Z’ in Eq. (A.7) is 
given by —Q!jj(0) + A where ^/^(O) = 0.5 is the dominant reggeon intercept and 

* The only factors which are not determined in this way are the factor [1/3 + 
2/3(1 —Z)] in Eq. (A.6) and the factor 2/3 in Eq. (A.8). These factors have a trivial 
origin. Eor instance the factor in (A.6) results from the fact that the fragmentation 
of u and d quarks into a proton is the same at Z —> 0. However, at Z —^ 1 the 
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A = 2a'p\ ~ 0.5 is common to all fragmentation functions and results from the 
integration over transverse momentum. In the case of A production has a 

softer Z ^ 1 behaviour, with a power —q;^( 0) + A = 0.5, which explains the different 
shape of A and p (non-diffractive) rapidity distributions (see Figs. 1 and 6). Of 
course, there are some numerical uncertainties in the determination of these powers 
- in particular in the power of Z in the leading baryon fragmentation function. The 
value 1.5 used above has been chosen, within the theoretically allowed range, in 
order to obtain the best description of the pp data in the DPM formalism. 

The only remaining point is to define the functions Gqq and Gq in the r.h.s. of 
Eq. (A.2), in terms of the fragmentation functions defined above. This is done in the 
following way. Let us consider first the production of antibaryons and non-leading 
baryons. In this case we have 

Gj,, 2 (y) = Sft z+ D!j{Z.) , for !, > A 

GJ, 2 (!/) = Z--Dj(Z-) -D,\, 2 (Z+) . for!,<A . (A.9) 

Here Z^ and Z^ are the light cone momentum fractions defined as Z^ = exp(y®^^ — 
Vm’ax)^ = exp(-y®^^ - V^max)^ where = y - A and V^Fax is the maximal 
rapidity that a hadron h can have in the string (for fixed values of the momentum 
fractions Xi and X 2 of the string ends). Note that for a string with a very large in¬ 
fragmentation of a M-quark is three times larger since the three combinations u{ud), 
u{du) and u{uu) are possible. (The later will give a A++ which will often decay into 
a proton). On the contrary, for the fragmentation of a d quark the only possible 
combination is d{uu). Likewise, the fragmentation of a diquark uu into a leading 
neutron is only possible through diquark breaking and has a probability of 2/3 as 
compared to the corresponding one for a ud diquark. 
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variant mass Z- ~ 0 for y > A and therefore G^q 2 {y) Z^ Dqq{Zj^). Likewise, 

for y < A, ~ 0 and ~ Z_ D^{Z_). 

In the case of a leading baryon, Eq. (A.2) has to be modified. We have 


dy 


(y A, Sgtr) — Gqq^l(y) 


(A.2') 


for both y > A and y < A, with 


= a, Z+ . 


(.4.10) 


when the diquark is located in the positive rapidity hemisphere. In the opposite 
case one has to replace Z^ by Z-. Note that, when the string invariant mass is 
large, one has Z^ i{Z+) ~ 0 in the backward hemisphere and, therefore, the 
leading baryon distribution is concentrated on the diquark side - as it should be. A 
further modification is necessary in order to preserve baryon number conservation. 
Eqs. (A.l) (A.2’) do not satisfy this property - a main source of violation being 
the threshold factor. Actually, a non Monte Carlo formalism which leads to exact 
baryon conservation has not been found so far. However, a good numerical accuracy 
in this conservation law is achieved by dividing the r.h.s. of Eq. (A.l) by 

dxi / dX 2 p{xi) p{x 2 ) 0 {Sstr - Sthr) ■ (A. 

Jo 

This device to enforce baryon number conservation must be used only for leading 
baryon production. 

The normalization factors ah and dh have been determined in order to repro¬ 
duce the measured rapidity distributions of baryons and antibaryons. We have : 




ttp = 0.82 , ttp = 0.036 , ttA = 0.17 , ay = 0.012 


(A12) 
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Note that the same normalization constant a is used for antibaryons and for non¬ 
leading baryons. 

Plugging the above expressions for the rapidity distributions of the individual 
strings in Eq. (2.3) we can compute the baryon and antibaryon rapidity distribu¬ 
tions in nucleus-nucleus collisions. The only missing ingredients are the fragmenta¬ 
tion functions of the strange sea diquarks us, ds and ss. For these functions we use 
the following approximate relations 


= dN = 

■^us,l -^ds,l -^us,l -^dsA ^udA ’ 


-^SS,l -^SS,1 -^UU,1 


(A.13) 


Although everything is now fully especified in order to compute baryon and an¬ 
tibaryon rapidity distributions, we give, as an illustration, the detailed form of the 
fragmentation of a qqg-Qs string into a leading A (which is, of course, identical to 
that of an A in a qqs-qs string). We have 

= ^ ^ pI\{x 2) e {Sstr - Sthr.) {A.U) 

where 

Gqq^liy) = ^2 + + 

KdAZ+) + . (A.15) 

Eq. (2.4) is obtained from (A.l) after dividing by the denominator in (A.11). Note 
that the p function of the sea diquark is taken to be the same as the corresponding 
one of a sea quark. In Eq. (2.5) we have explicitly written the contributions of 
the non-strange sea diquarks (terms proportional to ua) and that of the strange 



24 


sea diquarks (terms proportional to Up). The contribution of the non-strange sea 
diquarks (mm, dd, ud and du) is proportional to 4/(2 + S)‘^. That of the diquarks 
containing one strange quark (ms, ds, su and sd) is proportional to 4S/{2 + 5')^, 
and that of the diquark ss to S^/{2 + 5')^. For the fragmentation functions of the 
strange sea diquarks we have used Eqs. (A. 13). 

Note that for the baryon production in strings containing sea diquarks we have 
used the same prescription as for valence diquarks, namely we have divided by 
the denominator D (see Eqs. (A.11) and (A.14)) - although, in this case, this is 
not required by baryon number conservation since baryons are produced by pairs. 
Without dividing by D the production of A and A from sea diquark strings would 
have been substantially smaller. Actually, we regard the hadronization of the sea 
diquark strings as the main source of uncertainty in this approach. 


Turning next to the production of p and p, we have to replace in Eq. (A. 14), 


Gqq,l by 


Glqpiy) = 


4 


(2 + S) 


2 


|dL,i(z+) + 


(A.16) 


In this case, the contribution coming from the fragmentation of strings containing 
strange sea diquarks turns out to be very small and has been neglected. 


In order to compute the final state interaction of secondaries (pions and nu¬ 
cleons) introduced in Section 3, we also need the rapidity distribution of h~. In 
this case the contribution of the q-q strings is not negligeable and has the effect of 
narrowing the rapidity distribution. Eor our present purpose, however, it is enough 
to consider the two valence strings, provided the fragmentation function of valence 
quark and diquark are chosen in such a way as to reproduce the NN h~ dis¬ 
tribution. As expected from the absence of q-q strings, the diquark fragmentation 
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functions have to be softer than expected from Regge arguments. (Of course this 
drawback disappears when multistrings are introduced [32]). We take 

GlAy)=a^-(i-Zf ■, Gl^^a^Al-Zf (A17) 


where 

Z = |2yU7,sinh(y - A)/^^| , 

and /Utt is the transverse pion mass. The normalization, chosen to reproduce the 
rapidity plateau height in NN collisions, is a^- = 0.59. For the quark fragmentation 
functions we take the ones introduced in Refs. [33], namely 

Gf (!/) = -^(l + Z) F(Z) 

of {y) = ^{1 - Z) F{Z) (A18) 


where 


F{Z) = 


1.3(1 - + 0.05 

1 - 0.5Z 


Actually, in computing the h rapidity distribution in A A collisions, we have 
used Eq. (2.2), i.e. the expression without sea diquarks. More precisely, we take 

dN^fAT 


dy 


= nA (+f+”(!/) + Jv£7/(!/)) + 


(n - ua) 


2 + S 


(y) + NF a- (y) 


-P T 
>ls -Is 


(A.19) 


In Eq. (A. 19) we only consider the contribution of the non-strange sea strings. Of 
course some pions are also produced in strange sea ones. However their contribution 
is small and roughly compensates the loss of pions resulting from the introduction 
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of sea diquark strings (see Section 2), where pion production is smaller than in Qs-Qs 
strings. 

Finally, to compute the kaon multiplicities (Table 1), we have used the exper¬ 
imental values in NN collisions [6] and have assumed 11] that the fragmentation 
s (s ^ K~) is equal to u ^ {d ^ 7r“). This is, of course, justified only 

at the first string break-up. However, due to the small momentum fraction carried 
by the sea quark the probability of further break-ups in the rapidity hemisphere of 
the sea quark is negligeably small. 
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Figure Captions 

Fig. 1 Feynman x distributions of A and A produced in pp collisions at 200 GeV/c. 
The compilation of experimental data, at energies ranging from 100 GeV/c 
to ^/s = 63 GeV, is from Ref. 20. Error bars in the A data are only shown 
for X = 0.2 and 0.95. 

Fig. 2 Rapidity distribution of A in central SS collisions at 200 GeV/c per nucleon 
(dashed line) compared to the data of Ref. 4. The dashed-dotted line is 
the corresponding result without final state interaction. The full line is the 
prediction for central Pb-Pb at 160 GeV/c per nucleon, scaled down by a 
factor 5. 

Fig. 3 Rapidity distribution of A in central SS collisions at 200 GeV/c per nucleon 
(full and dashed lines) compared to the data of Ref. 4. The dashed-dotted 
line is the corresponding result without final state interaction. The dashed 
(full) line is obtained by adding the final state interaction, Eq. (3.3), with 
the proton spectrum corresponding to the dashed (full) lines in Eig. 7. The 
upper lines are the predictions for central Pb-Pb collisions at 160 GeV/c 
per nucleon scaled down by a factor 5. The dashed (full) line is computed 
by adding the final state interaction with the proton spectrum computed in 
DPM (shifted by Ay = 0.5). 
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Fig. 4 The rapidity distribution of negative hadrons {h ) in central SS collisions 
(full line) and in N-N collisions (dashed line), the latter scaled up by a factor 
32, compared to the data of Ref. [34]. 

Fig. 5 Rapidity distribution of negative hadrons {h~) in central Pb-Pb collisions at 
160 GeV/c per nucleon. The data are from Ref. [35]. 

Fig. 6 Feynman x distributions of proton and antiprotons in pp collisions at 200 
GeV, compared with the data of Ref. [36] at 175 GeV/c. Diffraction is not 
included in the model calculation. 

Fig. 7 Rapidity distribution of proton minus antiproton in peripheral (full line) and 
central (full and dashed lines) SS collisions compared with data from Ref. 
[17]. The theoretical curve for the peripheral case corresponds to N + N ^ 
{p — p) + X, normalized to the data. The dashed line in SS is the result 
of DPM. The full line in SS is obtained by adding some extra stopping to 
the proton (with Ay = 0.5) without changing the proton average multiplicity 
(see main text). 
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Table Caption 

Average multiplicities of secondaries in NN, central SS collisions at 200 GeV/c 
per nucleon and central Pb-Pb at 160 GeV/c per nucleon, compared with avail¬ 
able data. For nucleus-nucleus collisions the multiplicities of A, and Kq in 
DPM are also given before the final state interactions tiN —> iFA and tiN KK 
are taken into account. The number of protons in Pb-Pb (followed by an as- 
terix), corresponds to the combination {p n)/2. The values of multiplicities 
within (without) a parenthesis correspond to the dashed (full) lines in Fig. 3. The 


nucleon multiplicities after final state interaction are obtained from the relation 


DPM+f.s.i _ DPM _ ( DPM+f.s.i _ DP_m\ 
>(p) -^P{P) 'v A(A) ^A(A) ) 


1.6/2 (see main text). 



















